Anomalously strong pinning of the filling factor v = 2 in epitaxial graphene 
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We explore the robust quantization of the Hall resistance in epitaxial graphene grown on Si- 
terminated SiC. Uniquely to this system, the dominance of quantum over classical capacitance in 
the charge transfer between the substrate and graphene is such that Landau levels (in particular, 
the one at exactly zero energy) remain completely filled over an extraordinarily broad range of 
magnetic fields. One important implication of this pinning of the filling factor is that the system 
can sustain a very high nondissipative current. This makes epitaxial graphene ideally suited for 
quantum resistance metrology, and we have achieved a precision of 3 parts in 10 10 in the Hall 
resistance quantization measurements. 

PACS numbers: 73.43.-f,72.80.Vp,06.20.F- 



X 



The quantum Hall effect (QHE) is one of the key fun- 
damental phenomena in solid-state physics It was 
observed in two-dimensional electron systems in semicon- 
ductor materials and, since recently, in graphene: both in 
exfoliated an d epitaxial [H-Q devices. A direct high- 
accuracy comparison of the conventional QHE in semi- 
conductors with that observed in graphene constitutes a 
test of the universality of this effect. The affirmative re- 
sult would strongly support the pending redefinition of 
the SI units based on the Planck constant h and the elec- 
tron charge e [Io| and provide an international resistance 
standard based upon quantum physics (llj . 

Graphene is believed to offer an excellent platform for 
QHE metrology due to the large energy separation be- 
tween Landau levels (LL) resulting from the Dirac-type 
"massless" electrons specific for its band structure [12{ . 
The Hall resistance quantization with an accuracy of 3 
parts in 10 9 has already been established Q in Hall-bar 
devices manufactured from epitaxial graphene grown on 
Si-terminated face of SiC (SiC/G). However, for graphene 
to be practically employed as an embodiment of a quan- 
tum resistance standard, it needs to satisfy further strin- 
gent requirements (llj . in particular with respect to ro- 
bustness over a range of temperature, magnetic field 
and measurement current. A high measurement current, 
which a device can sustain at a given temperature with- 
out dissipation, is particularly important for precision 
metrology as it defines the maximum attainable signal- 
to-noise ratio. 

The extent of the QHE plateaux in conventional 2D 
electron systems is, usually, set by disorder and temper- 
ature. Disorder pins the Fermi energy in the mobility gap 
of the 2D system, which suppresses dissipativc transport 
at low temperatures over a finite range of magnetic fields 
around the values corresponding to exactly filled LLs. 



These values can be calculated from the carrier density 
n s determined from the low-field Hall resistivity measure- 
ments and coincide with the maximum non-dissipativc 
current, the breakdown current. Thus, the breakdown 
current in conventional two-dimensional scmicondutors 
peaks very close to the field values where the filling fac- 
tor v is an even integer [ll|. Though less studied ex- 



perimentally, the behaviour of the breakdown current on 
the plateaux for the exfoliated graphene, including the 
v = 2 plateau corresponding to the topologically pro- 
tected N = LL, looks quite similar (l3T |. 

In this Brief Report we explore the robustness of the 
Hall resistance quantization in SiC/G. Unlike the QHE 
in conventional 2D systems, where the carrier density is 
independent of magnetic field, here specifically to SiC/G, 
we find that the carrier density in graphene varies with 
magnetic field due to the charge transfer between surface 
donor states in SiC and graphene. Most importantly, we 
find magnetic field intervals of several Tesla, where the 
carrier density in graphene increases linearly with the 
magnetic field, resulting in the pinning of v = 2 state 
with electrons at the the chemical potential occupying 
SiC surface donor states half-way between the N = 
and N = 1 LLs in graphene. Interestingly, at magnetic 
fields above the v = 2 filling factor pinning interval, the 
carrier density saturates at a value up to 30% higher 
than the zero-field carrier density. The pinned filling fac- 
tor manifests itself in a continuously increasing break- 
down current toward the upper magnetic field end of the 
v = 2 state far beyond the nominal value of -B„=2 cal- 
culated from the zero-field carrier density. Facilitated by 
the high breakdown current in excess of 500 fj,A at 14 T 
we have achieved a precision of 3 parts in 10 10 in the Hall 
resistance quantization measurements. 

The anomalous pinning of v = AN + 2 filling factors in 
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SiC/G is determined by the dominance of the quantum 
capacitance, c q , [l4j over the classical capacitance per 
unit area, c c , in the charge transfer between graphene 
and surface-donor states of SiC/G: c q <C c c , where c q = 
e 2 7 e , c c = l/(4nd) and "f e is the density of states of 
electrons at the Fermi level. The latter reside in the 'dead 
layer' of carbon atoms, just underneath graphene floT - l^ . 
This layer is characterised by a 6\/3 x 6-\/3 supercell of 
the reconstructed surface of sublimated SiC. Missing or 
substituted carbon atoms in various positions of such a 
huge supercell in the dead layer create localised surface 
states with a broad distribution of energies within the 
bandgap of SiC (« 2.4 eV). 

It appears that the density of such defects is higher 
in material grown at low temperatures (1200 — 1600°C) 
resulting in graphene doped to a large electron density, 
n a ~ 10 13 cm' 2 , which is difficult to change [2l|. On the 
other hand, growth at higher temperatures, T rj 2000°C, 
and in a highly pressurised atmosphere of Ar seems to 
improve the integrity of the reconstructed 'dead' layer, 
leading to a lower density of donors on the surface and, 
therefore, producing graphene with a much lower initial 
doping [71. 122} . 

The quantum capacitance of a two-dimensional elec- 
tron system is the result of a low compressibility of the 
electron liquid determined by the peaks in 7 e . For elec- 
trons in high- mobility GaAs/AlGaAs heterostructures in 
magnetic field, the quantum capacitance manifests it- 
self in weak magneto-oscillations of the electron density 
[23L I24I ] due to the suppressed density of states inside 
the inter-Landau level gaps. A similarly weak effect has 
been observed in graphene exfoliated onto n-Si/Si02 sub- 
strate 25( , where the influence of a larger (than in usual 
semiconductors) inter-LL gaps is hindered by a strong 
charging effect determined by a relatively large thickness 
of Si02 layer. For epitaxial graphene on SiC, due to 
the short distance, d « 0.3 — 0.4 nm, between the dead 
layer hosting the donors and graphene, the effect of quan- 
tum capacitance is much stronger, and the oscillations 
of electron density take the form of the robust pinning 
of the electron filling factor. A similar behaviour was 
observed in STM spectroscopy of turbostratic graphite, 
where charge is transferred between the top graphene 
layer and the underlying bulk layers 
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FIG. 1. (Color online) Schematic band-structure for graphene 
on SiC in zero field (a) and in quantising field (b); graphical 
solution for carrier density as a function of magnetic field, 
n a (B), of the charge-transfer model given by Eq.[T](black line) 
together with lines of constant filling factor (red/gray lines) 
and n s (B,N) (green/light lines) for n g = 5.4 x 10 11 cm~ 2 
(c) and n s = 8.1 x 10 11 cm -2 (d). The vertical dashed line 
indicates the maximum field of 14 T in our setup and the blue 
dot indicates v = 2 calculated from n s (0). 



transfer in SiC/G is illustrated in the sketches in Fig. [T] 
for B = (a) and quantising magnetic fields (b). The 
transfer can be described using the charge balance equa- 
tion EDl 



n s , of this charge density is transferred to graphene, and 
an amount, n„, (controlled by the gate voltage) - to the 
polymer gate [22(. 

Graphical solutions for the charge transfer problem for 
two values of n g are shown in Fig. [TJ (c) and H{d) for 
a broad range of magnetic fields. For graphene within 
26(. The charge interval III (visible only in the case of the higher n g ) 
the Fermi energy coincides with the partially filled zero- 
energy LL, ep = 0, which determines the carrier den- 



l+e 2 j/c c 



and can be up to 30% higher 



j[A - 4ire 2 d(n s + n g ) 



(1) 



The left hand side of this equation accounts for the deple- 
tion of the surface donor states, where A is the difference 
between the work function of undoped graphene and the 
work function of electrons in the surface donors in SiC, 
ep is the Fermi energy of electrons in graphene, and 7 is 
the density of donor states in the dead layer. An amount, 



sity n 

than the zero- field density n s (0) in the same device [21 
This regime of fixed electron density is terminated at 
the low field end, at Bm — hn QCl /2e : where the N = 
LL is completely occupied by electrons with the density 
rioo. Note that for the n g presented here, Bm > 14 T 
- the maximum field in our setup. Similarly, for mag- 
netic field interval I, the Fermi level Ef = HvV2/\b coin- 
cides with the partially filled N = 1 LL (As = y/h/eB), 



and, for this interval, n\ 



,(B,1) with n s (B,N) 
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n oo — lh \+ei^l c ~ ■ The interval I is limited by the field 
values for which the N = 1 LL in the electron gas 
with the density n l s is emptied at the higher field end, 
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g 1+e 2 7 / c j ■ In magnetic field interval II the chemical 
potential in the system lies inside the gap between N = 
and N = 1 LL in graphene. As a result, over this entire 
interval the N = LL in graphene is full and N = 1 is 
empty, so that the filling factor in graphene is fixed at 
the value v = 2, and the carrier density increases linearly 
with the magnetic field, n s = 2eB/h, due to the charge 
transfer from SiC surface. 

According to Eq. [1] lowering the carrier density using 
an electrostatic gate is equivalent to effectively reducing 
the work function difference between graphene and donor 
states by n g (l/j + e 2 /c c ), which shifts the range of the 
magnetic fields where pinning of the v = 2 state takes 
place. For instance, reducing the zero-field carrier density 
from n a = 6.7 x 10 11 cm" 2 [Fig. [TJc)] to n s = 4.6 x 
10 11 cm" 2 [Fig. [T] (d)] moves interval II from 11.5 T < 
B n < 21.6 T down to 7.7 T < B n < 15.9 T, almost 
entirely within the experimental range. 

In order to verify the predictions of the theory regard- 
ing the pinning of the v = 2 filling factor and its implica- 
tions for the resistance metrology, we studied the QHE in 
a polymer-gated epitaxial graphene sample with Hall bar 
geometry of width W = 35 pm and length L = 160 pm. 
Graphene was grown at 2000° C and 1 atm Ar gas pres- 
sure on the Si-terminated face of a semi-insulating 4H- 
SiC(OOOl) substrate. The as-grown sample had the zero- 
field carrier density n s = 1.1 x 10 12 cm -2 . Graphene 
was encapsulated in a polymer bilayer, a spacer polymer 
followed by an active polymer able to generate acceptor 
levels under UV light. At room temperature electrons 
diffuse from graphene through the spacer polymer layer 
and fill the acceptor levels in the top polymer layer. Such 
a photo-chemical gate allowed non-volatile control over 
the charge carrier density in graphene. More fabrication 
details can be found elsewhere [71 122[. 

Figurc[2ja) shows magneto transport measurements on 
the encapsulated sample tuned to a zero-field carrier den- 
sity of n s = 6.7 x 10 11 cm" 2 corresponding to the case 
in fig. He)]. From the carrier density we estimate that 
the magnetic field B v= 2 needed for exact filling factor 
v = 2 in this device is 13.8 T. A well-quantized Hall 
plateau in p xy can be seen at v = ±2 for both magnetic 
field directions which is more than 5 T wide, whereas the 
longitudinal resistivity, p xx , drops to zero signifying a 
non-dissipativc state. In addition, a less precisely quan- 
tized plateau is present at v = ±6, for which p xx remains 
finite. 

Accurate quantum Hall resistance measurements re- 
quire that the longitudinal voltage remains zero (in prac- 
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FIG. 2. (Color online) (a) Transverse (pxy) and longitudinal 
(pxx) resistivity measurement. The horizontal lines indicate 
the exact quantum Hall resistivity values for filling factors 
v — ±2 and ±6. (b) Determination of the breakdown cur- 
rent, Ic, for three different measurement configurations ex- 
plained in legend, (c) High-precision measurement of p xy and 
Pxx as a function of magnetic field. Ap xy /pxy is defined as 
(pxy(B) - p xy (l4:T)) / p xy (14T) and Pxy(B) is measured rela- 
tive to a 100 Q standard resistor previously calibrated against 
a GaAs quantum Hall sample 0|. All error bars are la. 



tice, below the noise level of the nanovolt meter) to ensure 
the device is in the non-dissipative state, which can be 
violated by the breakdown of the QHE at high current. 
Figure [U[b) shows the determination of the breakdown 
current I c at B = 14 T on the v = 2 plateau. Here 
we define I c as the source-drain current, I sc i, at which 
V xx > 10 nV. We find for three different combinations of 
source-drain current contacts that the breakdown current 
for this value of n s is approximately 50 pA (note that I s d 
in a practical quantum Hall to 100 51 resistance measure- 
ment is w 25 pA [13]). The contact resistance, deter- 
mined via a three-terminal measurement in the nondissi- 
pative state, is smaller than 1.5 O. 

Figure [2jc) shows a precision measurement of p xy and 
Pxx for different magnetic fields along the v = 2 plateau. 
Note that this plateau appears much shorter in the mag- 
netic field range than that shown in Fig. [3] (a) because 
of the 200 times larger measurement current used in 
precision measurements. From this figure we determine 



4 




600 



400 o 
< 

200 i-h u 



15 



FIG. 3. (Color online) (a) Experimental p xx (black line) and 
p xy (red/gray line) together with the measured break-down 
current, I c (blue squares), (b) Hopping temperature, T* as 
a function of magnetic field. Inset: ln(a xx T) versus T _1//2 at 
13 T. Red/gray line is linear fit for 100 > T > 5 K giving 
T* « 12000 K. 



that the mean of Ap xy /p xy is —0.06 ± 0.3 x 10~ 9 for the 
data between 11.75 and 14.0 T, while at the same time 
p xx < 1 mVt. This result represents an order of mag- 
nitude improvement of QHE precision measurements in 
graphcnc, as compared to the earlier record [?|]. Not only 
is QHE accurate, but it is also extremely robust in this 
epitaxial graphene device, easily meeting the stringent 
criteria for accurate quantum Hall resistance measure- 
ments normally applied to semiconductor systems. 



Using the polymer gating method |22j, we further re- 
duce the zero-field electron density n s in graphene to cor- 
respond to the solution of the charge transfer problem in 
Fig. [ljd), i.e. down to 4.6 x 10 11 cm -2 as evidenced 
by magnetotransport measurements in Fig. (3ja). On the 
v = 2 quantum Hall resistance plateau we measure the 
breakdown current / c , defined above, as a function of 
the magnetic field. Unlike the conventional QHE ma- 
terials 11| , the breakdown current in Fig. E^a) continu- 



ously increases from zero to almost 500 fiA far beyond 
B u =2 ~ 9.5 T calculated from the zero-field carrier den- 
sity. This is a direct consequence of the exchange of carri- 
ers between graphene and the donors in the 'dead' layer, 
which keeps the N = LL completely filled well past 

B v= 2- 

The magnetic field range where the Fermi energy in 
SiC/G lies half-way between the N = and 1 LLs de- 
termines the activation energy h^/T/2v / Xb ~ 1000 K 
for the dissipative transport. For such a high activa- 
tion energy, the low-temperature dissipative transport is 
most likely to proceed through the variable range hop- 
ping (VRH) between surface donors in SiC involving vir- 
tual occupancy of the LL states in graphene to which 



they are weakly coupled. Indeed, as shown in the inset 
of Fig.[3jb), the temperature dependence of the conduc- 
tivity a xx measured at B = 13 T obeys an exp(— ^JT* /T) 
dependence typical of the VRH mechanism. The T* val- 
ues determined from the measurements at different mag- 
netic fields are plotted in the main panel of fig-E^b). The 
breakdown current rising with field to very large values 
[Fig. Eta)] corresponds to T* reaching extremely large 
values in excess of 10 4 K - at least an order of magni- 
tude larger than that observedjn GaAs [28[ and more 
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recently in exfoliated graphene 

In conclusion, we have studied the robust Hall resis- 
tance quantization in a large epitaxial graphcnc sample 
grown on SiC. We have observed the pinning of the v = 2 
state which is consistent with our picture of magnetic 
field dependent charge transfer between the SiC surface 
and graphene layer. Together with the large break-down 
current this makes graphcnc on SiC the ideal system for 
high-precision resistance metrology. 
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